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Abstract. In the recent paper (Casselman, 2001) I described how a number of ideas due 
to Fokko du Cloux and myself could be incorporated into a reasonably efficient program to 
carry out multiplication in arbitrary Coxeter groups. At the end of that paper I discussed 
how this algorithm could be used to build the reflection table of minimal roots, which could 
in turn form the basis of a much more efficient multiplication algorithm. In this paper, 
following a suggestion of Bob Howlett, I explain how results due to Brigitte Brink can be 
used to construct the minimal root reflection table directly and more efficiently. 

Introduction 

A Coxeter system is a pair (W, S) where W is a group with a set of generators S and 
relations 

(st) m ^ = 1 

for pairs s, t in S. These groups play an important role in mathematics far from visible in 
this simple definition. Much work has been done on them, particularly on finite and affine 
Weyl groups, but many phenomena involving Coxeter groups remain unexplained, and it is 
likely that computer explorations will be even more significant in the future than they have 
been so far. These groups become extremely complex as the size of S grows, however, and 
computer programs dealing with them must be extremely efficient to be useful. Conventional 
symbolic algebra packages normally fail to handle the difficulties satisfactorily. This paper 
is the second of a series in which I describe what may be even in the long term the most 
efficient algorithms to do basic computations in arbitrary Coxeter groups, at least on serial 
machines. These algorithms depend strongly on mathematical results of Brigitte Brink, 
Robert Howlett, and Fokko du Cloux. 

The first problem one encounters is how to multiply two elements of the group, which is by 
no means a simple matter. For familiar examples such as S n there are many satisfactory 
solutions, and for the Weyl groups of Kac-Moody algebras one can use the representation of 
the group on the root lattice to reduce many problems to integer arithmetic. But arbitrary 
groups, even sometimes the extraordinary finite groups H3 and H±, are computationally 
more demanding. Even with Weyl groups integer overflow may occur when working with 
the root lattice. In addition, solving computational problems in this domain leads to math- 
ematically interesting questions as well. 

The best combinatorial solutions — certainly, those of greatest theoretical interest — seem to 
be those in which elements of the group are represented in terms of products of elements 
of S. More explicitly, impose an ordering on S. Suppose w to be equal to the product 
S1S2 ■ ■ - s n as a product of generators in S, and for each i let Wi be the partial product 
si . . . Si. This expression for w is called its InverseShortLex normal form if for each i the 
element Si is least in S such that t(wiSi) < £(wi). It is the shortest representation of w 
least in lexicographical order when read in reverse. The multiplication problem can now be 
posed: given the normal form of w, what is that of sw ? The first general solution was implicit 
in work published by Jacques Tits in the mid 1960's. It had the virtue of being entirely 
combinatorial in nature, working only with the normal form of w, but it was extremely 
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inefficient, requiring both time and storage space roughly exponentially proportional, or 
worse, to the length n of the normal form of w. 

A second solution was proposed by Fokko du Cloux around 1990, but explicitly described 
only for finite Coxeter groups. It was this that I worked out in more detail in (Casselman, 
2001). It is mildly recursive in nature, of unknown theoretical complexity but in practice not 
unreasonably slow. It also has the tremendous virtue of requiring little machine memory. 
If one wants to be able to carry out just a few multiplications, it is likely that this method 
cannot be improved on. But if one wants to carry out a great many multiplications for a given 
Coxeter group, then it is probably best to build an auxiliary data structure that will make 
multiplication much simpler. A structure perfectly suited to this purpose is a table describing 
how the generators of S act on the minimal roots of (Brink and Howlett, 1993). As will 
be recalled later on, it can be used in a multiplication algorithm of complexity essentially 
linear in the length of w, allowing multiplication to take place without backtracking in a 
left-to-right scan of the normal form of w. 

In the paper (Casselman, 2001) I proposed that du Cloux's algorithm could be used as a 
kind of bootstrap to build this reflection table, but as I mentioned at the end of that paper 
it is not difficult to find groups where the number of minimal roots is rather small, hence 
the table itself also fairly small, but the task of building it in this way extremely tedious. 
This is unsatisfactory. 

After I explained this method to Robert Howlett, he and I experimented with other tools, 
and after a short while he suggested that the results of Brigitte Brink's thesis seemed to be 
applicable. In this thesis, almost all of which has been summarized in her published articles, 
she describes rather explicitly what the set of minimal roots looks like for any Coxeter group. 
However, it was not at all apparent at first glance, at least to me, that what she does is 
computationally practical. Howlett showed that my first impression was false by producing 
a program in the symbolic algebra package Magma that builds the minimal root reflection 
table of an arbitrary Coxeter group. This implementation was much slower than subsequent 
programs written by Fokko du Cloux and myself, but already dealt easily with a previously 
difficult group described in my paper. It might be worthwhile to say that the new algorithm 
depends on several crucial results of Brink's thesis, but not on her more or less explicit lists, 
interesting though they may be. 

The principal drawback of the new algorithm as the basis of multiplication is that it claims 
an amount of machine memory proportional to the number of minimal roots. The number of 
minimal roots can grow quite fast with the size of S, so this is not a negligible objection. It 
thus exemplifies the usual programming trade-off between time and space. Familar Coxeter 
groups are deceptive in this regard — for classical series of finite and affine groups the number 
of minimal roots is roughly \S\ 2 , whereas one of the first really difficult groups Howlett's 
program dealt with was an exotic Coxeter group with l^l = 22 and several hundred thousand 
minimal roots. (Howlett's Magma program took 69 minutes to handle it, he tells me, whereas 
a program in C written by Fokko du Cloux took 10 seconds on a comparable machine.) Of 
course, such groups are so complicated that there is not much exploration one can expect 
to do with them anyway, and in practice the number of minimal roots is only a relatively 
minor impediment. 

In the rest of this paper I will first recall the basic properties of minimal roots; explain 
how they are to be used in multiplication; summarize results of Brigitte Brink's thesis in 
the manner most useful for my purposes; and finally explain how to apply these results in 
practice to construct the minimal root reflection table. I should say that what I am going 
to come up with is somewhat different from Howlett's program, and also quite different 
from a related program of Fokko du Cloux. It uses a few more advanced features of Brink's 
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thesis than Howlett's program. It is less sophisticated than du Cloux's program, but has 
the advantage of being possibly clearer in conception. 

There is not much originality in this paper. Once Howlett had suggested using the results 
of Brink's thesis, the path to a practical program was probably almost determined. On the 
other hand, this path did not appear obvious to me until after a great deal of experiment, 
and I feel it will be useful to place here a record of the outcome. I should also say that the 
order of exposition and its occasional geometric emphasis are new. I will be satisfied if this 
exposition awakens interest in the work of Brink, Howlett, and du Cloux. This whole area 
of research is one of great charm, but it has apparently and unfortunately awakened little 
interest at large. This is likely due to the fact that it lies somewhere in the great no man's 
land between pure mathematics and practical computation — terra incognita on most maps, 
even in the 21st century. I must also record that both my programs and this exposition 
have benefited enormously from conversations with Fokko du Cloux. As I learned long ago 
to my own chagrin, competing directly in programming with him can be an embarrassing 
experience, but working alongside him is, on the contrary, at once pleasant and educational. 

1. Minimal roots 

The algorithms to be described are intimately related to geometric realizations of Coxeter 
groups. Following Brink and Howlett, I shall work here exclusively with the standard 
realization. To each element s of S is associated an element a s of a basis A of a real vector 
space V. An inner product is defined on V according to the formula 

a s »a t = - cos(ir/m Stt ) = -(C 2m +C 2 "m)/ 2 ( m = m s,t, = e 2m / 2m ) . 

In particular, a s »a s = 1, a s *a t = —1 whenever s and t generate an infinite group, and 
— 1 < a s »at < 1/2 whenever they generate a finite non-abclian group. A representation of 
W on V is determined by mapping elements of S to reflections: 

sv = v — 2 (v • a s )a s . 

Recall that the Coxeter graph has as nodes the elements of S, and an edge linking s and 
t if m Sl t, which is called the degree of the edge, is more than 2. In other words, the edges 
link s and t with a s • a t 7^ 0. In case a • ^ 0, I write a ~ 0. A simple link in the graph is 
one with degree 3, a multiple link one of higher degree. The formula for reflection amounts 
to the specification that if A = ^il then under reflection by s a all the coordinates of A 
remain constant except that indexed by a, which changes to 

(sA) a = -A a + ^[-2a./3] Xp . 

In calculations, it is frequently best to work with 2 a • rather than the dot-product itself, 
since it is at once simpler and more efficient to do so. 

This representation is faithful. Let 

C = {v I v • a s > 0} . 

The faces of this simplicial cone are the closures of the open simplices 

C T = {veV\(a,v)=Q for a e T, (a, v) > for a <£ T} 

where TCS. Thus C itself is C$ , and a face of codimension one is some C a . The union of 
the domains wC is a convex conical region called the Tits cone C of the realization. The 
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group W acts discretely on C, and the closure C of the open simplicial cone is a fundamental 
domain for this action. Each face of wC is the transform of a unique Ct; it is said to be 
labeled by T. 

A transform wa s of an element of A is called a root of this realization. The root hyperplanc 
A = is fixed by the root reflection s\ where s wa — wSaW^ 1 , the conjugate of an element 
in S. For any root A, define 



A root A is called positive if A > on C, negative if A < on C. Every root is either positive 
or negative, which means that no root hyperplane ever intersects C. A root is positive if 
and only if it can be expressed A = XLeA ^ a0i w ^ tn ^ a — ®- 

1.1. Proposition. Suppose that A and fi are positive roots. The following are equivalent: 
(a) A • jU = cos(ir k/m s ,t) with < k < m s ,t for some s, t in S; 



(c) the group generated by the root reflections s\ and s M is finite; 

(d) the hyperplanes A = and \i = intersect in the interior of C. 

proof. This is implicit in (Vinbcrg, 1971). That (a) implies (b) is trivial; that (b) implies 
(c) is a consequence of the simplest part of his discussion of pairs of reflections. If G is any 
finite subgroup of W then the G-orbit of any point in the interior of C will be finite, and 
since C is convex the centre of mass of the orbit will be in the interior of C as well. This 
shows that (c) implies (d). Suppose L = { A = 0} n = 0} to contain a point in the interior 
of C. The whole of C is tiled by simplices wCt- Points in the open cones wC are fixed only 
by 1 in W, and those in some wC a are fixed only by 1 and a single reflection. Therefore 
the linear space L must contain some wCt with \T\ = 2. This means that s and t can be 
simultaneously conjugated into Wt, and shows that (d) implies (a). □ 

1.2. Corollary, ff A and fj, are distinct positive roots such that |A»/x| > 1, then cither Ca>o 
contains C^>o orC^o contains C\>q. 

proof. If the intersection of A = and \i = with the interior of C is null, then the sign 
of n on the intersection of A = with that interior will be constant. Suppose it is negative. 
I claim that the region C M >o is contained in the region C\>o- If not, suppose that P is a 
point of C with (fi, P) > and (A, P) < 0. If Pq is any point in the interior of C, then the 
open line segment (Pq, P) is contained in the interior of C, it is also contained in the region 
C^^o, and it must intersect the hyperplanc A = 0. At that point /x will be non-negative, a 
contradiction. If the sign is positive, switch A and n.n 

Following Brink and Howlctt, I say that the positive root A dominates the positive root 
jj, if the region Ca>o contains the region Cp>o. Another way of putting this is to require 
that A > in every chamber wC on which fi > 0. Loosely speaking, a minimal root is 
one which dominates only itself. It is an important and remarkable fact about arbitrary 
Coxeter groups, first proved by Brink and Howlett, that the number of minimal roots is 
always finite. The significance of this result does not appear in familiar cases, since all roots 
are minimal when W is finite, and the minimal roots for affine Coxeter groups are rather 
simple in nature. 




{veC\ (X,v) =0} 
{veC\ (X,v) > 0} 
{veC\ (X,v) < 0} . 



(b) |A./i|<l; 
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Figure 1. The domains A > for the seven 
minimal roots of a Coxeter group with ma- 
trix 

"13 4" 

3 1 3 

4 3 1 

(Exhibited in a slice through the Tits cone.) 

In the rest of this section I recall from (Brink and Howlett, 1993) some basic facts about 
minimal roots. Proofs and statements will often be the same as theirs, and I'll refer to them 
when this happens. But in some places what I take from them is only implicit in their 
paper, and in other places I have given both statements and proofs a geometrical flavour, 
somewhat as in (Casselman, 1995). 

1.3. Lemma. (B & H 2.2(h)) Suppose X and /i to be positive roots with A dominating [i. If 
w/i is positive, then wX is also positive, and dominates w/i. 

proof. Under the hypothesis, C w \>o contains C w ^>q. If wo. is positive, the second contains 
C, and therefore wX > on C as well.n 

It is well known that s in S permutes the complement of a s in the set of all positive roots, 
or in other words that if s is in S and X a positive root, then either A = ot s and sX < 0, or 
sX > 0. For the minimal roots, there is a similar range of options: 

1.4. Proposition. If X is a minimal root and s in S then exactly one of these three options 
holds: 

• A = a s and sX < 0; 

• sX is again minimal; 

• sX dominates a s . 
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proof. If A ^ a > 0, then s a X > 0. If it is not minimal, say it dominates (3 > with 
s a X If [5 7^ a, then s a f3 > 0, and by the previous Lemma A dominates s a /3. Since A is 
minimal, this can only happen if A = s a /3, a contradiction. □ 

The minimal root reflection table is a matrix p of size \S\ x N, where N is the number of 
minimal roots. The entries are either minimal roots, or virtual minimal roots I arbitrarily 
label as G and ©. I set p(s,X) = p if sX = p is minimal, p(s, A) = if A = a s , and 
p(s, A) = ® if sX is a non-minimal positive root. This table is one of the fundamental data 
structures of a Coxeter group. 

The depth (5(A) of a positive root A is the length of the smallest w such that tt; _1 A < 0. 
Equivalcntly, it is one less than the length of the shortest gallery C = Co, C\, . . . , C n 
with C n in the region A < 0. Recall that a gallery is a chain of chambers with successive 
chambers sharing a common wall. The depth is also the smallest n with 

A = s 1 s 2 ■ ■ ■ s„_ia 

for some a € A. The elements of A itself, for example, have depth 1. A partial order is 
thus induced on the set of all roots: A ^ p if p = wX where S(p) = £(w) + S(X). This order 
gives rise to the root graph, whose nodes are positive roots, with edges A — > sX whenever 
A -< sX. 

1.5. Lemma. If A > and s«A dominates a then A -< s a X. 

proof. If (d) is a gallery from C to the region where s a X < 0, there will be some pair Cj, 
Ci+\ sharing a wall on the hyperplane a = 0. The reflected gallery will be shorter. Hence 
S(X) < S(s a X).n 

The minimal roots lie at the bottom of the root graph, in the sense that: 

1.6. Proposition. (B & H 2.2(iii)) If X < p and p is minimal, so is X. 

proof. This reduces to the case p — s a X. If p is minimal but p ^ a, then either s a p = X is 
minimal, or it dominates a. But by the preceding Lemma, in the second case 5(p) < S(s a p), 
a contradiction. □ 

1.7. Lemma. (B k H 2.3) If X> and A • a > then s a X -< X. 

proof. Suppose ui _1 A < with £(w) = S(X). If w~ 1 a < then we can write 

w = s a s 2 . ■ ■ s n 
A = 3 a 3 2 ■ ■ ■ s n -ia n 
s a X — s 2 ■ ■ ■ s„_ia„ 
S(s a X) = n — 1 < S(X) = n . 

So we may assume that w~ 1 X < but w~ 1 a > 0. Let w have the reduced expression 
w = si . . . s n , and let u = ws n — si . . . s n -\. Thus wa n < 0. I claim that u~ 1 s a X < 0, 
which implies that 5(s a X) = n — 1. 

We can calculate that 

w~ l s a X = w~ 1 X — 2 (A • a)w~ 1 a 
u _1 s Q A = it _1 A — 2 (A • a)s n w~ 1 a . 

From the first equation, we see that w _1 A will be negative under the assumption that 
X*a > 0. Since the depth of A is n, p = u _1 A will be positive. But then p > 0, s n p < 
implies that p = a n . So we have 

u~ 1 s ct A = a n — 2 (A • a)s n w~ 1 a . 
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The root 

a n — 2 (A • a)s n w~ 1 a 

is either positive or negative. In the second term, the positive root w~ 1 a cannot be a n , 
since then wa n < 0. Therefore s n w~ 1 a is a positive root also not a n , and the second term 
is a negative root not equal to a multiple of — a n . The sum has to be negative. □ 

1.8. Corollary. (B & H 1.7) Whenever A is a positive root and s in S there are three 
possibilities: (a) sX — > A, (b) sX = A, or (c) A — ► sX, depending on whether X»a s > 0, 
A • a s — 0, or X • a s < 0, respectively. 

1.9. Corollary. (B & H 1.8) If X — J2 ^ A* — S C« a ^ ej3 ^ A 1 ^- 

These results will soon be used to outline an algorithm to construct the minimal root 
reflection table. Processing roots will take place in order of increasing depth, applying this 
criterion: if A ^ a is a minimal root, then s a X is no longer a minimal root if and only if 
X»a s < — 1. Details follow, after I explain in the next section how to use the minimal root 
reflection table in multiplication. 

2. Minimal roots and multiplication 

Suppose that w has the normal form w — S\ . . . s n . Then the normal form of sw will be 
obtained from that of w by insertion or deletion: 

NF(sw) = si . . . SitSi+i ...s n 

or 

NF(sw) = si . . . Si_is i+ i ...s n . 

How do we find where the insertion or deletion occurs, and which of the two it is? If 
insertion, what is inserted? 

The algorithm explained in (Cassclman, 2001) handles these problems in geometric terms. 
I define the InverseShortLex tree to be the directed graph whose nodes are elements of 
W, with a link from x to y if NF(y) — NF(x)»s. Its root is the identity element. Each 
y 7^ 1 in W is the target of exactly one link, labeled by the s in S least with ys < y. Paths 
in this tree starting at 1 match the normal forms of elements of W. 
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This tree gives rise to a geometrical figure in which we put a link from xC to yC if there is 
an edge in the tree from x to y. The multiplication problem can now be formulated: given 
an InverseShortLex path from C to wC, how can we End the InverseShortLex path 
from C to swC? As explained in (Casselman, 2001), the InverseShortLex figure is very 
close to being symmetric with reflection to the reflection s; only links near the hyperplane 
a s = are not preserved upon reflection. The chambers principally affected are what I call 
exchange sites of the InverseShortLex tree — wC is an exchange site if the chamber wC 
has a face on the hyperplane a s — 0, labeled by t in S which is less than the labelling of the 
InverseShortLex link entering wC. 
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Figure 2. Some of the exchange nodes for reflec- 
tion by (1) in A^. Links modified under reflection 
are colored. The fundamental chamber is at the 
bottom. 



The relationship between this notion and multiplication is simple. Suppose NF(w) = 
S\. . . s n and for each i let Wi be the partial product s\ . . . s,. Suppose that w k C is the last 
exchange site in the gallery C — Co, • ■ ■ , C n — w n C ' , and suppose that t labels the wall on 
a = 0. Then cither the gallery crosses a = there or it does not. In the first case, t = Sk+i 
and 

NF(sw) = s\ . . . SfcSfc+2 ...s n . 

In the second, 

NF(sw) — si . . . s k ts k+1 ...s n . 

The question is now: How do we recognize InverseShortLex exchange sites? This is where 
roots and, even more pleasantly, minimal roots come in. 

2.1. Lemma. Suppose w = Si . . . s n . The chamber wC is an InverseShortLex exchange 
site if and only if — /? with sp < s n . 

proof. The wall of wC labeled by j3 is embedded in the hyperplane w(3 = Q.n 

This suggests the following algorithm to calculate NF(sw): Read the string s± . . . s n from 
left to right. As we read, keep track of an element t in 5, an index k, and a root A. The index 
k is that of the last exchange, the element t is the label of the last exchange, A = w^ 1 A. To 
start with k = —1, t = s, A = ct s . As we read s,, A changes to s^A. If A = (3 with sp < Sj, 
t changes to sp and k to i. When we are through scanning, we have the data determining 
the last exchange. 

So far, minimal roots have not come into play. But in fact, in calculating the reflected root 
sA we only need to apply the minimal root reflection table. If A is minimal and sA = 0, 
then the chambers Wi-\Ci-\ and C; = Wi-\SiCi lie on opposite sides of the hyperplane 
a = 0, so we know we are looking at a deletion. If sA is minimal and sA = ® then we know 
the gallery has crossed a hyperplane A = which separates wC from a = 0. The exchange 
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we last recorded will be the last we ever record, because the gallery cannot recross that 
hyperplane to reach a = 0. 

We now know how to use the minimal root reflection table to multiply, and it remains to 
tell how to construct it. 

3. A rough guide to constructing the minimal roots 

A straightforward procedure to list all minimal roots, along with a description of how the 
elements of S act on them, is not difficult to sketch. 

The data determining a minimal root A are an index, a coordinate vector, and a list of 
reflections sX. There are at least two reasonable choices of coordinates. The coordinates I 
use in this note are those in the expression of a root as a linear combination of the roots in 
A: 



The other reasonable choice is the array (X*a), which du Cloux has discovered to be a 
choice somewhat more efficient if also more difficult to work with in programming. Indices 
are assigned in the order in which the minimal roots are first calculated; the roots in A 
are assigned the first |5| integers. Furthermore, in the preliminary version of the algorithm 
we shall maintain a look-up table of some kind, enabling us to find a root's data given its 
coordinates. 

A minimal root is to be dealt with in two stages: (1) defining and (2) finishing. In the 

definition, a root is calculated, assigned an index, entered in the look-up table, and put on 
a waiting list. Items are taken off this list in the order in which they are put on — the list 
is, in programmers' terminology, a FIFO list (first in, first out) or queue. 

A root A is finished when it is taken off the queue and all of the reflections sX not already 
known are found. At this moment, all roots with depth less than that of A have been finished, 
so that in particular we have recorded all reflections /j, = sX -< A. We must calculate the 
sA not yet determined; the problem is to figure out whether sX = A, sA is a minimal root, 
or sX = ®. Of course we can tell by inspection whether sX = A. Otherwise we know that 
A -< p = sX. There are three possibilities: (1) We have already defined /i, in which case it 
will be registered in our look-up table. We just add to its record that sX = p and sp = A. 
(2) The root p is not minimal. We can tell this when we calculate ^t, since in this case 
X*a s < —1. We set sX = ©. (3) Otherwise — 1 < A»a s < 0, and the root p is a new 
minimal root not yet defined. We do so by assigning it an index, putting it in the look-up 
table, recording the reflection sX = ^t, and putting p in the process queue. Then we go on 
to the next root in the queue, as long as there is something there. 

There are a number of problems that occur in trying to implement this is in practice. 
The first is that maintaining the look-up table is somewhat cumbersome. As du Cloux 
first observed, one of the results from Brink's thesis allows us to dispense entirely it. The 
following is a endsatz of the main Theorem of (Brink, 1995). 

3.1. Lemma. Suppose s =/= t in S. Assume the group W s j they generate to be Rnite. 
Suppose A to be any root whose W s . t -orbit doesn't intersect A. Then W S)t acts freely on 
this orbit. 

How can this be applied in our algorithm to find all minimal roots? Suppose we are in the 
process of defining a minimal root u,, which we may assume not to be dihedral. We calculate 
it as fi = sX with A -< p. If £ is an element of S with tX ~< A, then W Syt will necessarily be 
finite. In the orbit of p under W s . t there will be a unique element v of least depth, equal to 




Minimal roots 



11 



w s t /i where w s .t is the unique longest element of the dihedral group W s j- We can calculate 
it by finding alternately tX, stX, etc. until the chain starts to ascend in depth. And the 
difference in depth between /i and v will be exactly m s ,t precisely when tfj, -< /z. The element 
t/j, will then be {tw s .t)v, which can be calculated from the reflection tables of elements of 
depth less than A. In other words, as soon as u. is defined we can calculate all the descents 
t[i -< [i — without doing any new calculations! 

In the algorithm above we can therefore eliminate reference to the look-up table. The 
procedure sketched above requires, however, that we know whether the reflection of a root 
has less depth or not, so we add to the data of a root A its descent set, the set of s with 
sX -< A. When we define a root, before we put it in the queue, we calculate all of its descents. 

Because the result above is applicable only to roots not in the orbit of W s ,t, a U the dihedral 
roots have to be handled somewhat specially — they should all be defined immediately after 
the roots in A, and those of depth greater than 2 have to be fed into the queue at the right 
moment. 

There is another simplification. The following was observed by Brink and Howlett, and used 
as an important part of their proof that the set of minimal roots is finite: 

3.2. Proposition. If sX = © and X -< n. then s/i — © as well. 

proof. For a, [3 in A, A any root 

s a X»/3 = (A- 2(X»a)a) »/3= (A«/3) - 2(X»a){a»/3) 

which means that for (3 ^ a the dot product A»/3 decreases under reflection by a, as long 
as A < s a X. Therefore once we have calculated that sX — © this remain true for all \i with 

A < /!.□ 

I call an s with sX — © a lock at s. This Proposition tells us that locks are inherited as we 
go up the root graph and these inherited locks, like descents, can be immediately recorded 
upon definition. 

Eliminating the look-up table is a great simplification, noting inherited locks a minor one. 
In simple cases — when none of the Coxeter matrix entries m s j are greater than 3 — there 
is probably no better algorithm than the one sketched above, but incorporating the simple 
modifications concerning descents and locks. In this case, all doubled dot products 2 (A • a) 
will be integers. When A ~< s a X this doubled dot-product will be a negative integer, and if 
it is anything other than —1 the root sX will no longer be minimal. What could be simpler? 

If there are entries m s .t > 3 this still works, at least in principle, although there are a few 
annoyances that arise. All dot products in the range (—1,0) will be of the form cos (Trk/m s j) 
for k < TO^t/2. So in order to tell whether sX dominates a s we just have to compare the 
dot product X»a s to numbers in this finite set. There is one potentially nasty problem, 
however. All we can assert a priori about the coordinates of A is that they will be in the 
real cyclotomic field generated by cos(27r/M), where M is the lowest common denominator 
of the numbers 2m Sj t. The amount of work involved would seem to grow disturbingly as M 
grows. It is Brink's thesis that tells us that this apparent difficulty is spurious. In the next 
section I follow her in examining more carefully the coefficients of minimal roots. 
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4. Root coefficients 

We want to know something about what possible values can occur for the coefficients of 
roots and, especially, minimal roots. Roots are produced by a sequence of reflections from 
the basic roots in A. As is the case in many algorithmic processes, although we understand 
each single step — here, reflection — quite well, the overall development is not so clear. 

We begin with a discussion of dihedral roots. 

Suppose S to possess two elements s and t, corresponding to roots a and (3. Let m — m s , t - 
What are all the positive roots of the system? 

Let z = ( 2m = e" /m = e 2 "/ 2m . Then 2 a . (3 = 2 cos(tt/to) = (z + z" 1 ) and 

sA = A — 2(A • a)a 
2 (pa + q/3) • a = 2p - (z + z~ x )q 

s(pa + q/3) = [ — p+(z + z~ 1 )q]a + q/3 
tX = A- 2(A«/3)/3 
2 (pa + q(3)*{3 = 2q-(z + z'^p 

t(pa + q(3) =pa+ [ - q + (z + z~ l )p] (3 . 

Suppose m = oo. We get in succession (in tabular form) as transforms of a 

coefficient of a coefficient of (3 

a 1 

to 1 2 

sta 3 2 

(st) n a 2n + l 2n 

t(st) n a 2n + l 2n + 2 

and similarly for the transforms of (3. Of these, only a and (3 are minimal roots. 

If m = 2, the positive roots are a and (3. 

Now suppose 2 < m < oo. Then we get 

coefficient of a coefficient of (3 

a 1 

to 1 z + z^ 1 

sta z 2 + z~ 2 + 1 z + z^ 1 

(st) n a C n+1 C n 

t(st) n a C n+1 C n+2 

where the C n are solutions of the difference equation C n+ 2 = — C n + (z + z _1 )C„ + i with 
initial values Co = 0, C\ = 1, giving 

Cn = Z —^- = + z- 3 + • • • + z-(- 3 > + = sin( ; n / m) . 

z — z 1 sin(7r/m) 

Note that C n > C n -i as long as the real part of z n is positive, which happens as long as 
n < m/2. There are m roots in all for a dihedral system, which agrees with this remark. 

For m — 3 we get the table 
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For to = 4: 



a 
P 

ta — s(3 



a 
P 

ta = sta 
sp = ts/3 



For m = 5, with 9 = (1 + \/5)/2: 



a 

/3 
ta 
sp 

sta — ts/3 

For to = 6: 



coefficient of a 
1 

1 



coefficient of a 
1 

1 

V2 



coefficient of a 
1 




ta - 

s/3 



a 

P 
sta 
ts/3 



sta = tsta 
tsP = stsp 



1 



coefficient of a 
1 

1 

V3 
2 

V3 



coefficient of P 

1 
1 



coefficient of P 

1 

V2 
1 



coefficient of (3 

1 
1 
9 

e 



coefficient of P 

1 

V3 
1 

V3 
2 



The table for to > 6 looks somewhat similar to the last, except that it is longer. All 
coefficients beyond the first 6 rows are greater than 2. 

We can use these results, following Brink's thesis, to prove: 

4.1. Proposition. The coefficients of the positive roots are polynomials in the constants 

/-n n 

^2 (n < m/2) 

^2m ~ (2m 

with coefficients among the natural numbers N. 

proof. By induction on the depth of a root. If the depth is 1 there is no problem. Otherwise 

A 



suppose 



wa 



si . . . s n _ia t 



with a in A, (5(A) = n. Let s = s n . Choose y =/= 1 maximal in W s j with w — xy. Then 
xa s > 0, xa t > 0. We know from the discussion of dihedral roots what the coefficients of 



ya t = pa s + qa t 



are like, but then 



wa t — pxa s + qxa t ■ 
We can apply an induction hyothesis to xa s and xa t since £(x) < £(w).n 
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4.2. Corollary. The smallest positive root coefficient is 1. Any root coefficient lying between 
1 and 2 equals 2 cos(ir/m s ,t) for some s, t. 

PROOF. If 

sin(7m/m) 

• , I \ (! < n < TO /2) 
sm(7r/77j) 

lies between 1 and 2 then n = 2, and the least possible value is \/2. Any root coefficient 
will be a sum of terms 

C = ci . . . c fc 

where fc > 1 and each Cj will be either a positive integer or one of the constants in the 
Proposition, with each q > 1. The smallest possible value for one of these constants is 
2cos(7r/4) = \/2, so if this product is less than 2 then k = 1, and is an integer only if the 
product itself, of course, is l.n 

4.3. Corollary. The smallest possible coefficient of a positive root, other than 1, is y/2. 



5. Composition and decomposition 



One of the important results of Brink's thesis is that the apparent difficulties concerning 
mixed cyclotomy in the determination of the minimal roots are in fact spurious. In order 
to see that it is true, we have to look carefully at how roots are calculated. 

A root ^2 X a u may be considered as a function on the Coxeter graph: s i— > A Qa . The 
support of a root A — i s support as a function — the set of s where X as ^ 0. 

Brink's nice idea is to track the computation of roots by looking at their support. All roots 
are constructed by applying some w to a basic root. The support of a root A is extended 
by applying a reflection s which is not already in its support. If s is not connected to that 
support in the Coxeter graph, then sX = A. This process gives a new root only when s is 
linked to the support of A. Therefore: 

5.1. Proposition. The support of any root is connected. 

Suppose we apply a reflection s to a positive root A whose support docs not contain s. Let 
a = a s . Only the coefficient X a changes, and to 



If A is a minimal root with A Qs = 0, then sX will also be a minimal root if and only if 



[-2a'0\\f,<2. 



But Xp > 1 and also — 2 a* (3 > 1. The condition for minimality cannot hold here unless 
there is just one f3 = at in the sum. Furthermore, in that case the above condition amounts 
to the condition 

[-2a»/3]Xf3 < 2 . 

The first term is either 1 or > y/2, as is the coefficient A^. All in all: 

5.2. Proposition. Suppose X to be a minimal root, s an element of S not in the support O 
of X. If sX is again minimal then there must be just a single link from s to 6, say to t, and 
one of these holds: 
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(a) m Sj t = 3 and A Qt < 2; 

(b) 3 < m s ,t < oo and A Qt = 1. 

We can summarize this in a mnemonic diagram: 

• x- — ■ i — ► x- — x (x < 2) 

• 1=0 i — ► l=c (c = 2cos(7r/m)) 

5.3. Corollary. The support of a minimal root is a tree with containing no links of infinite 
degree. 

One of Brink's good ideas is to focus on the coefficients of a minimal root which are equal 
to the minimum possible value 1. Suppose A to be a minimal root with A Q = 1. We know 
that A has been built up from a chain of minimal root predecessors 

Ai -<...-< A„ = A . 

Some one of these, say \i = Afc, will be the first to have a non-zero coefficient at a, and that 
coefficient must be 1, since it can only grow. Unless fi = Ai = a, the root fj, must have been 
constructed by extension from a root Afc_i whose support did not contain s, and the earlier 
Proposition tells us that \i = sXk-i, where s is linked to the support of A^^i at a single 
node t by a simple link in the Coxeter graph, and \k-i,u t = ^k,a t — 1- We have A = wfi 
where w is a product of generators other than s. 

Since the support O of A is a tree, the complement of a in is the disjoint union of 
exactly as many components 0* as there are edges in the support of A linked to a, say I. 
Reflections corresponding to nodes in different components will commute. Hence we can 
write X — wi . . .wgfi where each Wi is in Wq' . All the Wi commute, so we can write the 
product in any order. Each of the Aj = Wia is itself a root with support in 0, = 9* U {a}, 
minimal because it will be a predecessor of A in the root graph. We have therefore proven: 

5.4. Proposition. Let X be a minimal root with X a = 1. Let be the support of A, the 0* 
the connected components of 0\{a}, 0j = 0* U {a}. The restriction Xi of X to 0^ is also 
a minimal root. 

This is called the decomposition of A. Conversely: 

5.5. Proposition. Let be a subset of S, a an element of 0, the 0* the components of 
0\{a}, 0j = 0* U {a}. If for each i the minimal root Xi has support on 0j with coefficient 
1 at a, then the vector X with support whose restriction to each 0j agrees with Xi is a 
minimal root. 

This is called composition. 

proof. More difficult than the previous result. It follows from this result of Brink's: 

5.6. Lemma. If X is any positive root with coefficient 1 at a, then X — wet, where w is a 
product of s in S not equal to s a . 

proof. Induction on the depth of A. Trivial if A = a. It remains to be shown that if 
<5(A) > 1 and A Q = 1 then we can find s such that sX -< X and sX also has coefficient 1 
at a. There exists at any rate some s with sX -< X. If this does not have coefficient 1 at 
a, then — as we have seen above — sa has coefficient at a; a has exactly one link to the 
support of sX, say to f3; and the coefficient of sX at j3 is equal to 1. By induction, sX = wj3 
where w is a product of t equal neither to s nor sp. But then 



A = swf3 = wsa .□ 
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These results have a natural generalization. The unit support of A is the set of a with 
X a = 1. Suppose A to be a minimal root, 6 its support, T C 6 its unit support. We know 
that 9 is a connected tree. The complement 0\T can be expressed as the disjoint union 
of its connected components 0*. For each of these, let 0j be its union with elements of T 
attached to it. The set ©i will then be connected, and no point of T will be in its interior. 
The set T will contain in addition some pairs of nodes of the Coxeter graph connected by 
links in it; let these be counted as additional sets Oj (but with empty interiors). Finally, 
there is one exceptional case where has just one element, which is ©i. 

5.7. Lemma. The set © is the union of the 6j, which overlap only in points ofT. 

proof. We may suppose © to have more than one element in it. Suppose given t in 6. If 
it is not in T, it belongs to a unique 0*. Otherwise, it will be in T. If it is connected in the 
Coxeter diagram to one of the 0*, it will be in 0j. Finally, we have the situation where the 
only elements of it is linked to are in T. Each of the pairs forming such a link will be one 
of the ©i.n 

If T is any subset of S, a subset is called T-connected if none of its interior points belong 
to T. Each of the 0, defined above arc T-connected, and are called its T-components. 
Following Brink, I call a minimal root indecomposable if its unit support is contained in 
the boundary of its support. 

Lemma. (Brink) (Decomposition) The restriction of a minimal root to any of the T- 
components of its support is again a minimal root. 

proof. This follows by induction from the earlier decomposition. □ 

5.8. Proposition. (Brink) (Composition) Let be any connected, simply connected subset 
of S, T an arbitrary subset of 0, 0, its T-components. If Aj is for each i a minimal root with 
support 0, and unit support Tn 0,, then the vector A with support whose restriction to 
0j is Aj will be a minimal root. 

This reduces the calculation of minimal roots to the calculation of indecomposable ones, 
together with a bit of accounting. Finding the reflection table turns out not to be much 
more difficult. 

How are indecomposable roots found? They are obtained by reflection from predecessors in 
the root graph. What can we say about this process? If A is an indecomposable root and 
ji = s a X -< A, what can we say about \xl Let be the support of A. There are several 
possibilities. (1) Suppose that a is a boundary point of 0, and that X a = 1. Then A is 
obtained from /x by extension from a unit node of fi. This must be a boundary point of the 
support of fi. This can happen only if fi has singleton support. (2) Suppose that a is in the 
boundary of but that A Q > 1. Then n, will also be indecomposable. (3) Suppose a lies in 
the interior of and fi a > 1. Then fi is also indecomposable. (4) Finally, suppose a lies in 
the interior of but \i a = 1. The support of fi is also 0, so it cannot be indecomposable. 
It will the composite of several \ii adjoined to each other at a. Thus: 

5.9. Proposition. An indecomposable root is constructed through a sequence made up 
of three types of processes: (1) extension, where fj, -< A = s^i, fi is also indecomposable, 
and the support of A is one node larger than the support of fi; (2) fusion, where [i is the 
composite of indecomposable roots in all with a in the unit support; (3) promotion, where 
fi is indecomposable and has the same support as A. 

We have already investigated extension. In the next section we look at fusion. 
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6. Reflecting at junctions 

Suppose to be a connected and simply connected subset of S, a € 0, the 0^ the connected 
components of 0\{s} joined with {s}. Suppose that A is a minimal root with support 
and A Qs = 1. Let Aj be the restriction of A to 0j, and assume the Aj are indecomposable. 
Under what circumstances is sX a minimal root? 

Each of the 0^ contains a unique node U linked to s. Let = ot ti . 

The inner product A • a is 

1- £>A \a'0i\ . 

Recall that Xp i > 1, |a«/3,| > 1/2. When the junction has more than 3 roots joining at it, 
this dot-product is < — 1, and sX will not be minimal. 

6.1. Proposition. If 0\{s} has more than 3 components, sX — 0. 
We may now assume there to be 2 or 3 links at a. 

6.2. Lemma. Suppose c — 2 cos(w/ m Si t)- The links c 1 and 1 = 1 cannot occur in any 

root. 

This introduces notation I hope to be self-explanatory. What I mean by this is that if A 
is a root then there do not exist neighbouring nodes s and t in the Coxeter graph with (1) 
m s ,t = 3, Xa s = 1, 1 < A Ql < 2; or (2) m s ,t > 3, X 0s = 1, A Ql = 1. 

proof. By induction. The only way to get c 1 is from c 1 l.n 

We now consider the possibilities for fusion. 

• Two links. 

The three a priori possibilities are 

X x ■ — ■ 1- — -X y 
X x ■ — ■ 1 = Xy 

X X = 1 = Xy 

6.3. Proposition. In the third case sX = ©. 

proof. The inner product is 1 — c x X x /2 — c y X y . We know that each c* is at least y/2, and 
by the previous Lemma that A* > \pl as well. Thus the dot product is —1 or less.n 

In the first case, the condition is 1 — Az/2 — A y /2 > —2, or X x + X y < 4. The case A* = c 
is excluded, so either A, = 1 or A, > 2. It is not possible for both A* to be > 2, so one at 
least must be 1, say X y = 1. 

We can do better: the inner product must be — cos(nk/m) for k < m/2. This gives us: 
m = 3. The inner product is 0, with 

• 1 1 1 i — > 1 1 1 

or —1/2 with 

• 2- — -1- — -1 i — > 2- — -2- — 1. 



m = 4. We must have 1/2 - A^/2 = -y/2/2 or A = 1 + V2/2. Impossible. 
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m = 5. 1 — A x = — c or 1 — c, A^ = 1 + c or A = c. The second is impossible, alas. So the 
only possibility is 

• (1 + c) — 1 1 !—» (1 + c) — (1 + c) 1. 

to > 6. Nothing possible. 
Now the remaining case 

X X 1 = Xy , 

Here 

1 - X x /2 - {c/2)X y > -1, A x + cA y < 4 . 
m = 4. Only possibility is 

• 1 1 = V2 i — > 1 2=V2. 

to = 5. Only possibility is 

• 1 l=c i — ► 1 (l + c)=c. 

Note that in all these cases we have a join of two primitive roots, one of which is just 1 1. 

So it is trivial to reconstruct the reflection even when this is all embedded in a generic root. 

• Three links. 

Assume the reflection is a primitive root. The conditions are X x + \ y + X z < 4 if all links in 
are simple, or X x + X y + c m X z < 4 if one has degree m > 3. More of degree to > 3 cannot 
occur by simple calculation. If to > 3 then c > \[2 so the second sum is at least 4, and 
cannot occur. 

As for the first, we cannot have any A* > 2, so each can only be 1 or c. But the link c 1 

cannot occur. 

6.4. Proposition. The only triple junction producing a minimal root be reflection at the 
junction is the case that all coefficients are 1. It reflects to a central node with value 2. 



i > 



1 1 I 1 

Figure 3. The only triple junction giving rise to 
fusion. 
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Looking back over all ways indecomposable minimal roots are constructed, we deduce: 

6.5. Proposition. No indecomposable minimal root has more than one link of degree m > 3 
in its support. 

This means that mixed cyclotomy does not occur. Brink's thesis at this point goes on to 
build a rather explicit description of all possible indecomposable roots whose support is a 
given Coxeter graph. It is probably a poor idea to try to use these lists in a computer 
program, since it is fairly simple to use basic properties of indecomposable roots to build 
these lists automatically. There is one consequence of what we have seen so far that is simple 
to formulate. It turns out that any root coefficient may be described in no more than three 
integers: each one is either of the form a + bc m with c m — 2 cos(7r/m) or larger than 2. With 
a bit of fiddling to handle the last case carefully, the three integers (a,b,m) are sufficient. 
This is a drastic version of the principle of 'no mixed cyclotomy'. Implementing this idea in 
an explicit program is straightforward enough, but the program to be sketched in the next 
section, which uses composition and decomposition in a more sophisticated fashion, seems 
to be more efficient. 

There is one final point to take into account. When an indecomposable root A has degree 
m > 6, it will have a unique link in its support of that degree. The discussion of extension 
above tells us that this link is in the support of the dihedral root first encountered in 
constructing A. Suppose this dihedral root has support at s and t, and that it started with 
s. The discussion of dihedral roots tells us that the coefficient of a s in A will be either 1 or 
greater than 2. In the first case, s has to be on the boundary of the support of A because of 
indecomposability. In the second, we know that it is on the boundary because no node with 
coefficient > 2 can be extended. In all cases, therefore, s is on the boundary of the support 
of A. I leave it as an exercise to show further that in the second case the coefficient of ct t 
is equal to c m — 2co&(tt /m s _ t ), and in both cases all other coefficients are of the form nc. m 
where n is a positive integer. (This, too, is an observation of Brink's thesis.) 

7. The program 

I am going to sketch here a program that lists all minimal roots and constructs the reflection 
table. Relying strongly on Brink's thesis, it works with two kinds of root, indecomposable 
and composite. The associated data structures are quite different. Both have indices, 
descent sets, and a reflection list. But an indecomposable root has in addition an actual array 
of coefficients. These coefficients will all lie in one of the real cyclotomic rings generated by 
C 2m + Ci2to> where to = m Sjt for some s, t. I call these rings elementary. An indecomposable 
root is also assigned its degree to, and when to > 6 its two special nodes, the ones spanning 
the unique link of degree to, are specified as well. A composite root has, instead, a list of 
its indecomposable components. 

New roots are found by applying a reflection to a root already defined. There are several 
basic ways in which to do this: (1) extension of an indecomposable root to produce a new 
indecomposable root with a larger support; (2) promotion of an indecomposable root to 
produce a new indecomposable root with the same support; (3) fusion at a junction of 
several indecomposable roots to produce a new indecomposable root; (4) composition of an 
arbitrary root with a dihedral root attached at one end to its unit support; (5) replacement 
of one or more of the components of a composite root by a single indecomposable root. 

It is only the first three that require actual arithmetic, and always in one of the elementary 
real cyclotomic rings. Composition is an extremely simple operation, while replacement 
amounts to reducing a composite reflection to one of the first three. That this can be done 
easily seems almost an accident depending on the limited number of ways fusion can take 
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place. I'll not give details of this below, but I'll sketch an example here. Suppose we want 
to calculate S3 A where 

A: 1 1 1=V2 V2 

and s = S3 is the third node from the left. The root A is the composite of three indecom- 
posable roots 

1- — -1, 1- — -1, 1=V2- — -V2 

There are two components in the star of s, the components of A containing it, and they do 
not exhaust A. Instead, they make up the smaller root 

p. 1 1=V2 V2 

whose reflections we shall have already calculated when A is taken off the queue. It is easy 
enough to locate /x, given the reflection tables already made, since it is equal to S2 applied 
to the indecomposable root 

v: 1=V2 V2 

We look up s<iv to get /j,, then look up sfi, which is the indecomposable root 

1 ■ — ■ 2 =V2- — -V2 . 

Finally, we replace the two components of [i by the single component sy, in the list of 
components of A. All this can be carried out, with minor modifications, for each of the cases 
of fusion. In the pseudo-program below the details will be grossly telescoped. 

Here is the program: 

First of all we define the basic roots, those in A. 

Then we define specially the dihedral roots, those whose support consists of exactly two 
nodes of the Coxeter graphs. When we are through with these we have finished the basic 
roots. We put the dihedral roots of depth 2 in the queue, and after we have defined all of 
those we put those of depth 3 in the queue. In this way, we preserve the essential property 
of the queue that roots of depth d are put on before those of depth d+1. 

As for dihedral roots of depth 4 or more, we can finish them as soon as they are constructed. 
Both of their coefficients are 2 or more, which means they cannot be extended, so the only 
reflections higher in the root graph are to dihedral roots with the same support. In other 
words, we never put them in the queue. Incidentally, these occur only for links with m > 6. 
One consequence is that in calculating reflections of roots taken off the queue, we have only 
to do arithmetic with algebraic integers of the form a + bc m , where m is the degree of the 
root, c m = 2cos( / n/m) 1 a and b are integers. 

At this point, all roots of depth 1 have been finished, all dihedral roots have been defined, 
and some have been finished. The queue contains only dihedral roots. While the queue is 
not empty, we remove roots from it to be finished. 

Let A be a root removed from the queue. All roots of less depth have been finished already, 
and in particular all descents sA -< A have been assigned. We run through all the s in 5, 
and for each one of them where sA has not been assigned, we have to calculate sA. Exactly 
what this involves depends on a number of things — the basic idea is to use composition and 
decomposition to avoid repeating actual arithmetic. 

So now suppose we are looking at s and A with sA not known. Let a = a s . There are many 
cases to deal with. 
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o The node s lies in the support of A. 
o It belongs to a unique component A^. 

o That component is A itself, which is therefore indecomposable. We must calculate 

the dot product A • a to see if sA is minimal or not, and if it is define it. 
o The component Aj is not all of A, which is a composite root. We have calculated s\i 
already, and sA is obtained by replacing Ai with s\i. 
o It belongs to several components. It is therefore a junction, and we are looking at 
a possible fusion. Let \i be the star of s, the composite of those components whose 
supports contain s. 

o The root /i is the same as A. We have to calculate sA explicitly. It may be A, ffi, or 

a new indecomposable root, 
o The root /i makes up only a part of A. We have calculated s/x, and obtain sA by 
replacing /i with sfj,. Exactly what we do to find /i and s/i depends on which case of 
fusion we are dealing with, as I mentioned in the example I looked at earlier. This 
can be decided by local information around s. 
o The node s does not lie in the support of A. 

o It is linked to more than one node in the support, or by a link of infinite degree. Here 
sA = ©. 

o It is linked by a single link of finite degree. 

o The link is to a unit node of A. Then sA is a new composite, 
o The link is not to a unit node. 

o The root A is indecomposable. We have to calculate explicitly, 
o The root A is composite. If fi is the component it is linked to, we have already 
calculated s[i. We replace /i by sfi in A. 
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